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The recent Belle and BaBar measurements of the branching ratio of — > t^v-j- indicate 
a significant deviation from the standard model prediction. We demonstrate that this mea- 
surement has a serious impact on models with minimal flavor violation involving a charged 
Higgs boson, ruling out a large portion of the currently allowed parameter space. In the 
constrained minimal supersymmetric standard model, this creates a tension between the 
measurements of — > t^v^ and the anomalous magnetic moment of the muon, unless 
tan /? is small, /i > 0, and Aq takes a large negative value. In fact, a very small region of the 
parameter space of this model, with small values of mo and mi/ 2, survives all the constraints 
at 95% C.L. It is remarkable that this specific region is still consistent with the lightest su- 
persymmetric particle as the dark-matter. Moreover, it predicts observable supersymmetric 
signals in the early runs of the LHC, even perhaps at 7 TeV. We also show that a consistent 
explanation for the deviation of the B^ — )■ t^Vt branching ratio from the standard model 
can be achieved in a nonuniversal Higgs-mass model, which could also predict early signals 
of supersymmetry at the LHC. 

PACS numbers: 12.60.Jv, 13.20.He, 14.80.Da 



I. INTRODUCTION 

Now that the CERN Large Hadron Collider (LHC) has commenced its long-awaited run and the 
first physics results have been analyzed and made public , there is an atmosphere of palpable 
suspense in the high energy physics community as to what physics results will come out as more and 
more data are collected and studied, and most importantly, as to whether these results will indicate 
new physics (NP) beyond the Standard Model (SM). The experimental programme is more or less 
clear: more statistics will be accumulated, and the results will be compared with the predictions 
of the SM. Deviations from the latter would imply some sort of NP, and one can refer to existing 
theoretical studies to indicate what kind of NP is indicated by the observed deviation. It is true 
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that theorists have not succeeded in providing an unequivocal prediction in this regard. This is 
because there exist several rival possibilities for NP, each with good arguments both for and against 
it. However, for several technical and aesthetic reasons, of which tractable ultraviolet behavior and 
the natural appearance of chiral fermions are perhaps the most important, supersymmetry (SUSY) 
has always been the pick of these NP models. At the dawn of the LHC era, it still remains the 
first option for any study of NP predictions. 

Elegant as SUSY may be as an abstract idea, it is well known that it presents a very different 
face when it comes to constructing realistic models at low energies. Any phenomenologically viable 
SUSY model must necessarily include a large number of soft SUSY-breaking parameters. A count 
of the number of free phenomenological parameters in the so-called minimal supersymmetric stan- 
dard model (MSSM) [3i] runs to over 100, including masses, coupling constants, and mixing angles 
for the large number of supersymmetric partners, or sparticles, in the model. This proliferation of 
parameters may be directly traced to the fact that the MSSM does not include a specific mechanism 
for the breaking of SUSY, and hence, the numerous SUSY-breaking parameters are essentially put 
in by hand. Although such a model can exist, at least in principle, a theory with a hundred odd 
free parameters is a phenomenologist's nightmare, since it leads to very few clear predictions at the 
empirical level. At the LHC, for example, this leads to a wide landscape of possible signals which 
would leave an experimentalist with hard data to compare with a bewildering variety of options 
J]. It is also difficult to believe that the breaking of SUSY is a sheer accident brought on by a 
proliferation of arbitrary nonzero parameters. One would rather argue that there is a definite mech- 
anism for SUSY breaking [5^, and when we know it, we will also know the parameters in question. 
Once again, however, theorists have failed to come up with an unambiguous mechanism for SUSY 
breaking - there exist quite a few different suggestions Q], beginning with minimal supergravity 
(mSUGRA) models, through gauge-mediated SUSY breaking (GMSB), anomaly mediated SUSY 
breaking (AMSB) and so on, each with a very different pattern for the parameters in question. 
Each of these models has different predictions for LHC signals, and hence, in effect, the chaotic 
situation within the subset of SUSY models becomes a cameo of the general NP scenario. 

The oldest, and perhaps the most restrictive, of these SUSY models where a specific mechanism 
for SUSY breaking is considered, is the so-called "constrained" MSSM (cMSSM), which is based 



on an underlying mSUGRA 



7|] scenarij^. In this scenario, supergravity is broken spontaneously in 



a so-called "hidden" sector consisting of fields which do not have strong or electroweak couplings 

^ This mechanism, by invoking supergravity, gives up the good ultraviolet behavior, unless, indeed, the supergravity 
model is equivalent to, or embedded in, a string theory. Aficionados of the cMSSM would, of course, argue that 
gravity must eventually be included anyway. 
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to the MSSM fields. However, gravity, which necessarily couples to all fields so long as they 
carry energy and momentum, acts as a mediator between the hidden sector and the MSSM sector, 
giving rise to the soft SUSY-breaking parameters. It is this circumstance that leads to a dramatic 
reduction in the number of parameters, since gravity is blind to all flavor and color quantum 
numbers, though it can sense the spin of a particle. As a result, the mSUGRA model has just 
five free parameters, viz. (i) a universal scalar mass tjiq, (ii) a universal fermion mass mi/2j (hi) a 
universal trilinear (scalar) coupling Aq, (iv) the ratio of vacuum expectation values of the two Higgs 
doublets, parametrized as tan /3, and (v) the Higgsino mixing parameter fj,. This universality of 
the masses and couplings is valid at the scale where the SUSY-breaking parameters are generated, 
which is usually identified with the scale of grand unification (GUT scale J^, i.e., above 10^^ GeV. 
While running down to low energies using the renormalization group (RG) equations, however, 
the various soft SUSY-breaking parameters evolve differently, and lead to a specific mass spectrum 
at the electroweak scale. In particular, one of the Higgs-mass-squared parameters is driven to a 
negative value, ensuring that the electroweak symmetry is spontaneously broken. The requirement 
that the electroweak symmetry be broken at precisely the right scale leads to a further constraint, 
which effectively fixes the magnitude of /x in terms of the other parameters, though its sign is still 
indeterminate. This version of the mSUGRA model, which depends on four parameters and a 
sign, viz. {mo, mi/2i ^Oi tan /3, sgn fi}, is called the cMSSM. Being more constrained, this model 
is also more predictive and hence is more readily testable. There exists, therefore, a vast amount 
of literature on this model, which has been studied with regard to (a) collider signals Q], (6) low- 
energy processes, such as decays of K,D, and B mesons P, [3], and (c) dark- matter constraints 
arising from the fact that the relic density of the lightest SUSY particle (LSP) can be identified 
with the dark-matter content of the Universe [Ul . In this article, therefore, we shall focus on this 
model, though a simple extension will also figure into our analysis. 

It is now common knowledge that null results from direct searches have pushed up the masses 
of sparticles into the regime of 100 GeV or above. However SUSY models can still make sub- 
stantial contributions to low-energy processes, particularly those which are mediated by weak 
interactions. Among these, flavor-changing neutral current (FCNC) processes, with the famed 
Glashow-Iliopoulos-Maiani (GIM) cancellation, constitute a favored ground to look for SUSY ef- 
fects (or any NP effects, for that matter). However, barring a few little hiccups, the SM rules 
supreme in the area of flavor physics, leaving very little room for NP theories, including SUSY and 



^ There exists a symbiotic relation between SUSY and grand unified theory (GUT) ideas, since SUSY solves the 
hierarchy problem in GUT , and a GUT is natural at the scale where SUSY breaking is generated. 
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the cMSSM. Year by year, as the measurements of the FCNC processes grow better and better, 
the lower bounds on masses of new particles (including sparticles) have been creeping further and 
further up in order to squeeze the NP contributions into the ever-narrowing band of experimental 
errors in these measurements. 

In this article, we consider one such recent low-energy experimental result, viz., the measure- 
ment of the branching ratio — )• t^l^t- It directly constrains all models with minimal flavor 
violation (MFV), viz., models where all flavor-changing transitions are entirely governed by the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix with no new phases beyond the CKM phase 5. We 
find constraints on general NP with MFV that involves a charged Higgs boson. As the cMSSM 
(and almost any viable SUSY model) belongs to this category, we apply these constraints to the 
cMSSM and find a rather dramatic impact on the parameter space of the model. It turns out when 
we combine the results of the measurement in question with other low-energy measurements, such 
as the anomalous magnetic moment of the muon, and radiative and leptonic B decays, most of 
the cMSSM parameter space is disfavored at the 95% confidence level (C.L.). What survives all 
the constraints is a small patch in the four-dimensional parameter space {mQ,mi/2, /3) of 

the model, for a positive sign of ^. This is very different from the kind of constraints derived from 
earlier, less restrictive measurements of — t- t^Ut, where wide areas of the cMSSM parameter 
space were allowed. 

As mentioned above, one of the attractive features of SUSY is that it provides a dark-matter 
candidate, viz., the lightest super symmetric particle (LSP). This carries a conserved quantum 
number (R parity) which forbids its deca}o. One can, therefore, study the evolution of the Universe 
in a SUSY model, and check whether the relic density of LSP's matches with the observed density 
of dark-matter as indicated by the cosmic microwave background radiation (CMBR) data [2]. 
Obviously, this matching will happen for only a small part of the parameter space of the model. It 
is encouraging that the dark-matter-allowed region in the cMSSM overlaps the small patch allowed 
by low-energy measurements quite substantially. We can say, therefore, that there exists a rather 
specific set of parameters which is simultaneously consistent with the low-energy data as well as 
with the hypothesis that LSP's form the dark- matter content of the Universe. With this specific 
set of parameters, we generate the mass spectrum of sparticles, and find reasonably unequivocal 
indications as to the kind of signals expected at the LHC. No detailed analysis is necessary at 
this stage, for the relevant signals have already been considered in comprehensive studies by the 
ATLAS and CMS Collaborations Q, [3]. Comparing their results with our parameter choice, we 
^ Once again, this is not written in stone, for R-parity violation can happen and has been extensively studied [l^ . 
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find that the 7 TeV run of the LHC may provide a weak indication of SUSY [l6^, which could be 
verified comprehensively even in the very early stages of the 14 TeV run. Going further, we may 
even say that if SUSY is indeed the correct NP option, then the LHC may eventually turn out to 
be the hoped-for SUSY factory, claimed in the literature fl?! . 

This article is organized as follows. In Sec. II we discuss the recent bounds on — t- t^i^-j- 
and how they affect MFV models. This is followed by Sec. Ill, where we discuss other low-energy 
measurements which constrain the cMSSM parameter space. The combined constraints are dis- 
played and discussed in Sec. IV, where we also discuss the possible LHC signals which could arise 
therefrom. Sec.V discusses the so-called nonuniversal Higgs-mass (NUHM) model, a variant of the 
cMSSM, and some of its consequences. A critical summary of our results forms the substance of 
the concluding Sec. VI. 



II. THE DECAY B+ r+Vr 

On purely theoretical grounds, the leptonic decay — t- r^Vr is a clean decay mode, since the 
final state consists only of leptons and hence the usually troublesome strong rescattering phases 
are absent. Indeed, in the SM, the branching ratio of B^ — )• t^Uj- is given by the tree-level formula 

BR(i?+ ^ r+.^sM = fl - 4V fl \Vu,?rB , (1) 

Svr V ^bJ 

where Gp is the Fermi constant, tb is the B^ lifetime, Jb = 192.8 it 9.9 MeV is the B^ decay 
constant, and mB,mr are the masses of i?^,r"'", respectively. Here 

iKfcl = (3.52 ±0.11) X 10"^ (2) 

is the relevant CKM matrix element, obtained through the combined fit [2^ to all the data 
excluding the B~^ — >• t^v-,- measurements. The SM prediction, including higher-order corrections, 
is 

BR(5+ ^ r+zy^)sM = (0.81 ± 0.15) x 10"^ . (3) 
As recently as 2008, the experimental average value of this parameter [2l| was 



BR(5+ ^ r+i/^)2008 = (1-41 ± 0.43) x 10"* , (4) 

which was just about consistent with Eq. Q at 1 standard deviation. At that time, it was shown 
22l | that corresponding constraints on the parameter space of the cMSSM (such as we discuss in 



this work) are rather minor. 
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Very recently (2010), however, new measurements of the branching ratio BR(i5"'" — ?■ t'^i't) from 
B factories have ch ang ed the experimental value quite significantly. The most recent experimental 



measurements are 



23|426| 



Babar (semileptonic tag) :BR(S+ T+Vr) = (1.70 ± 0.82) x 10~^ , 

(hadronic tag) :BR(S+ ^ r+Vr) = (1-80 ± 0.61) x 10"^ , 

Belle (semileptonic tag) :BR(S+ ^ t+v^) = (1.54 ± 0.48) x 10"^ , 

(hadronic tag) :BR(S+ ^ T+Vr) = (1-79 ± 0.71) x 10"^ . (5) 

These results are quite consistent with each other. Combining these measurements, one gets the 
world average j^] 

BR(S+ ^ r+i/^)exp = (1.68 ± 0.31) x 10"'^ . (6) 
Clearly, this measurement deviates significantly from the SM prediction given in Eq. ([3]). Defining 



-Rr^r to be 



28|, 



29| 



^ BR(i?+ ^ T+i^Oexp , . 

"^^-^ BR(i?+ ^ T+Z^OSM ' ^ ^ 

and using Eqs. ([3]) and ([6]), we get 

i?^^P = 2.07 ± 0.54 , (8) 

which indicates a ~ 2(T deviation. Deviations at this level frequently arise from statistical fiuctu- 
ations in small data samples, or from the use of ill-determined theoretical quantities, and often 
disappear when more data are analyzed, or when more rigorous calculations are performed. In this 
case, however, the mismatch may not disappear so easily. For neither are the current measurements 
of BR(5^ —7- T^Vr) based on a small statistics sample (see Refs. 23h26||). nor should one expect 



the SM prediction to change much, since the formula in Eq. ([T]) involves quantities that are already 
known pretty accurately. It appears, therefore, that it is sensible to at least explore the ability of 
NP beyond the SM to resolve the observed discrepancy. 

Following Ref. [28|], we characterize the NP models that could potentially explain this anomaly 
by a quantity -R^j^, defined as 

oNP ^ BR(^+ T+;^^)sM+NP 

~ BR(i?+ ^ T+Z.,)SM ' ^ ^ 

where the subscript SM+NP represents the net branching ratio in the NP scenario, including the 
SM contribution. The 95% C.L. allowed range for i?^,^ then works out to 

0.99 < < 3.14, (10) 
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which essentially means that NP models with positive contributions are favored by the data and 
those with negative contributions are quite strongly disfavored. 

There exist, of course, a wide variety of models of NP which could provide extra contributions 
to the branching ratio of — t- t^Vt- However, we focus only on the MFV models. For a large 
class of MFV models that involve a charged Higgs boson - such as two-Higgs doublet models, 
of which the cMSSM is an example - the branching ratio of — )• t'^Vt is given by [s^ 

BR{B^ ^ r+.^NP = - ^) ' \Vu,\' TB (l - tan^ ' (11) 

at the tree level, where M_|_ is the mass of the charged Higgs boson. Here NP stands specifically 
for the MFV model, but we retain the notation "NP" in the interests of simplicit}C. 

In the above formula \ Vub\ is the value of \Vub\ obtained in the context of the NP model, which 
in general will be different from \ Vuh\ obtained from the data in the context of the SM. In order 
to get rid of this uncertainty in the CKM parameter, we restrict ourselves to the determination of 
\Vub\ through only those measurements that are independent of NP. Such a fit is called the fit to 
the universal unitarity triangle (UUTfit) {si!], and it uses only 

• the measurements of |Ktf,/lc6| from semileptonic B decays, 

• the ratio of mass differences in the Bg and B^ systems: \AMs/ AM^l, and 

• the measurement of sin 2/3 from the time-dependent CP asymmetry in i?^ — t- J/il)K^*\ 
The UUTfit value of \Vub\ comes out as [l| 

iKbluuTfit = (3.50 ± 0.12) X 10-3 , (12) 

which is actually very close to the global fit in Eq. ([2]). Using this value, the SM prediction for the 
branching ratio of B~^ — t- t^Ut changes slightly from Eq. ([3j) and becomes 

BR(S+ ^ r+i^^)sM = (0.80 ± 0.15) x 10"^ . (13) 

Note that while the UUTfit Q is obtained using the lattice prediction fs = 200 ± 20 MeV {32] of 
the LQCD Collaboration, we use the more recent, averaged value from lattice simulations, fs = 



^ Note that our analysis for the MFV models in this section closely follows that of 19], with minor differences. 
Our constraints in the M+-tan /3 parameter space naturally are almost identical. However, we present the detailed 
analysis here for the sake of completeness and clarification of our procedure. 



192.8 lb 9.9 MeV [lo], which has a much smaller erroiO, for the calculation of BK{B~^ — )■ t+^'t-)sm- 
The 95% C.L. allowed range for R^^^ assumes the value 



0.99 < R^^^ < 3.19, 



(14) 



which forms the basis of all subsequent analyses in this article. Once \ Vub\ is chosen in this "model- 
independent" way, we can take \Vub\ = |^ub|uuTfit) and hence the theoretical MFV prediction for 



R^^ at the tree level becomes 
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If higher-order corrections are included then this ratio gets modified [33] to a form 
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NP 
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tan2/3 
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B 



(15) 



(16) 



1 + eo tan /3 

where eo encodes all the higher-order corrections, which, of course, will have some dependence on 
the free parameters of the MFV model. We take the range of eo to be —0.01 < eo < 0.01, as 
obtained in by a scan over reasonable values of the MFV model parameters. When a specific 
model, such as the cMSSM, is considered, eo can be calculated explicitly. 
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FIG. 1. (a) The dependence of R^J^ on the mass Af+ of the charged Higgs boson in MFV models for two 
values tan/3 = 10 and 50, and (b) the 95% C.L. constraints on the Af+-tan/3 plane. The vertically hatched 
regions in (a) correspond to higher order corrections varying between eo = —0.01 and +0.01, while the la 
(2cr) experimental measurements of R^J^ are shown by horizontal broken (solid) lines. The dark band in 
(b) corresponds to the LEP bound. The large, vertically hatched region in (b) is disallowed by the recent 
(2010) R^J' constraint, while the horizontally hatched region is disallowed by the 2008 data. 



^ Ideally, of course, the UUTfit needs to be performed again with the updated /a value. We have assumed that the 
updated fit will not significantly affect the |Kii,| value. 
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The impact of the experimental data on MFV models with a charged Higgs boson, as discussed 
above, can be clearly discerned from Fig. [U^a), where we plot the value of R^,^ as a function of 
the charged Higgs boson mass M^. As Eq. (jl4p indicates, such a model should tend to make 
ii^j^ greater than unity, and there is very little room for R^,^^ < 1. However, the negative sign 
on the right side of Eq. ()16|) indicates that unless the NP contribution is very large, the models in 
question have a tendency to diminish R^,^^ rather than augment its value. As a result, a model 
with a heavy charged Higgs boson cannot be considered as an explanation for the deviation of 
BR(i3"^ —7- T^Ur) from its SM value. Instead, if we do have such a model, we would expect rather 
strong constraints on its parameters, since the NP contribution must be squeezed into the small 
tolerance below unity, as given in Eq. (jl4p . Such a situation would naturally arise for large M_|_, 
when the NP contribution becomes negligible, and R^J^ — s- 1 for all tan /? - though it always stays 
less than unity. This corresponds to the rising part of the curves, towards the right end of Fig. [TJo). 
A glance at the figure will, however, leave no doubt that this limiting case is barely allowed at 2a 
for low tan /3 ( = 10), but disallowed for high tan /3 (= 50). We surmise, therefore, that for high 
values of M^, the — )• r^i/^ measurement favors low values of tan/3. 

For low values of M_|_, on the other hand, Eq. (jl6p tells us that it is possible for the NP 
contribution to be so large that it completely dominates the SM contribution, and in this limit, it 
is possible to have R^^^^ > 1. The explicit condition is 

1 + eo tan/3 ~ m| ' ^ ' 

However this can happen only for a very restricted set of (M_|_,tan/3) values, some of which are 
already constrained experimentally. For example, for tan/3 = 10, one can have R^,^ well inside 
the 2a range only if M4. is < 50 GeV, but such low M+ values are ruled out by the LEP direct 
searches, which give M+ > 79.3 GeV [35]. On the other hand, if tan/3 = 50, the same LEP data 
permit 140 GeV< < 220 GeV which can render R^J^ well inside the 2a range, as shown in 
Fig. [TJa). Thus, one may complement our earlier assertion by the statement that in the opposite 
limit, i.e. for low values of M+, the — t- t^Vt measurement favors high values of tan/3. This is 
also the limit in which the NP models contribute positively in accounting for the deviation of the 
experimental data from the SM predictions. 

The two limits are made explicit in Fig. W(J>), which shows the 95% C.L. constraints on the 
M_|_-tan /3 plane. The dark band represents the LEP constraint > 79.3 GeV and the vertically- 
hatched region is disallowed by the B^ — t- t'^Vt measurement. This leaves only two small unshaded 
regions for low M_|_ and high M_|_, in accordance with the above discussion. Our result may 
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be contrasted with the constraints obtained from the 2008 data, which are shown by horizontal 
hatching, and constrain only a small region with < 400 GeV and somewhat high tan /3. 

One may say, therefore, that the recent measurement of the branching ratio for — ?• t^Ut has 
considerably improved the constraints on MFV models with charged Higgs bosons. As the cMSSM 
belongs to this category, we should expect correspondingly severe constraints on the corresponding 
parameter space when we compare its predictions with this new experimental result. 



III. OTHER CONSTRAINTS 



When we consider an all-encompassing model like the cMSSM, with far-flung implications in 
almost all areas of electroweak physics, the constraints arising from — t- t^Ut cannot be consid- 
ered in isolation, but must be combined with other bounds - some of which are equally restrictive, 
at least at the 2a level. These constraints can be classified as follows. 

1. Theoretical constraints - arising from requirements of internal consistency of the model. In 
particular, if the Higgs-mass parameter which should be driven negative by RG running 
remains positive, we cannot explain electroweak symmetry breaking (EWSB) in this model. 
There is also a substantial region where the model predicts that the charged stau ti is the 
LSP, and is therefore precluded by the absence of a large relic density of charged particles. 



2. Collider bounds - arising from the nondiscovery in direct searches [35|] at the CERN LEP 
and Fermilab Tevatron of predicted particles, most notably the light Higgs boson and 
the lighter chargino xt- While the chargino couplings are large enough for the experimental 
bounds to practically saturate the kinematic reach of these machines, the light Higgs boson 
mass in SUSY models is generally sensitive to higher-order corrections, where there is a 
theoretical uncertainty of around 3-4 GeV at the next-to-next-to-leading order (NNLO) and 
higher |3fi]. To take care of this, we consider a softer lower bound of 111 GeV, rather than 
the kinematic bound of 114.4 GeV usually applied to the SM Higgs boson. 

3. Indirect bounds-arising from measurements of low-energy processes where new particles and 
interactions in NP models can also contribute. In the context of the cMSSM, the most 
important of these are the measurements of (a) the anomalous magnetic moment of the 
muon, (6) the rate of the radiative decay B^ — >• Xgj, and (c) the BR for the leptonic decay 
Bs — )■ n^n^. Here we have assumed that the NP is of the MFV kind, and that it survives 
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the measurements other than those exphcitly mentioned above. In particular, the large Bg- 
Bg mixing phase, or the Acp{B — )• Ktt) measurements, cannot be explained by any MFV 
models, and we assume that these anomalies will disappear with more data or with better 
theoretical calculations. 

Of the above, the theoretical and direct search constraints may be considered firm constraints, 
as they are unlikely to be changed by inclusion of further types of NP along with the cMSSM or 
whatever model is being studied. On the other hand, constraints from indirect measurements are 
not so robust, as they can easily change if some new effect is postulated. Before we proceed to 
apply these constraints to the cMSSM parameter space, therefore, a brief discussion of the actual 
measurements used in our analysis is called for. This forms the remaining part of this section. 

• The anomalous magnetic moment of the muon, = {g — 2)/2.' This is one of the most 
compelling indicators of NP and it is well known as a major constraint for NP theories such 
as supersymmetry or extra dimensions. The latest measured value [32] for is 

a^fP = (11659208.0 ± 6.3) x 10"^° . (18) 



In the recent past, the SM prediction j3ll] for a^j has undergone numerous vicissitudes with 
respect to the experimental data, occasionally being consistent with it and occasionally 
deviating at the level of 2c7-3cr. Much of the difficulty in making this prediction accurate 
lies in the fact that the experimental measurement is sensitive to two-loop corrections where 
some nonperturbative QCD effects due to the low mass scale are involved. The latter have 
to be obtained by fitting experimental data, whose errors then feed into the theoretical 
uncertainty. The most recent SM prediction is 38| 



a™ = (11659178.5 ± 6.1) x 10"^°. (19) 

The discrepancy between the SM and experiment is, therefore, 

Aa;f P = a;fP - af^ = (29.5 ± 8.8) x 10"^° . (20) 

This is at the somewhat high level of ~ 3.4(T, but is not normally considered a "smoking gun" 
signal for NP for reasons stated above. Nevertheless, in order to check if this discrepancy can 
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be explained with the cMSSM, we use a procedure [39(] that does not calculate the two-loop 



SUSY corrections, but includes them in the theoretical errors, to obtain a 95% C.L. range 

11.5 X 10~^° < Aa^^ < 47.5 x 10~^° , (21) 

where Aa^^ is the extra contribution due to NP. In our analysis, the NP in question will 
be the cMSSM, or a variant, but we choose, as in the previous section, to retain the label 
"NP". 

Obviously, in the cMSSM, the value of Aa^^ will depend on all the free parameters of the 
model. However, it is known that the sign of the cMSSM contribution is directly sensitive 
to the sign of the fi parameter [4^: for /i < 0, the cMSSM contribution is negative, while for 
^ > 0, a positive contribution is predicted by some regions of the cMSSM parameter space. 
Since the 95% C.L. range of AaJJ^ indicated in Eq. (f2T]) is entirely positive, it indicates 
that the sign ^ < is disallowed by the measurement of the muon anomalous magnetic 
moment, and even with ^ > 0, there are strong constraints on the remaining parameters of 
the cMSSM. 

The radiative decay — )• Xg^: In the SM, the BR of the radiative decay — )• has 
been calculated [4]| to NNLO in QCD to be 

BR(Brf X^7)sM = (3.15 ± 0.23) x 10"^ . (22) 

The current experimental average for the BR by the Heavy Flavor Averaging Group 



(HFAG) [42] is 

BR(Bd ^ X,7)exp = (3.55 ± 0.26) x 10"^ , (23) 

which is consistent with the SM prediction within 1 standard deviation, leaving very little 
room for NP contributions. As a result, this measurement has a tremendous impact on 
MFV models involving a charged Higgs boson ^ essentially pushing up the mass M+ 
to very large values. Of all such models, the constraints on SUSY models can be more 
relaxed because of large cancellations between the charged Higgs boson contributions and 
the chargino contributions which are the hallmark of SUSY models. Nevertheless, there do 
exist bounds on the cMSSM arising from the residual contribution, especially for light M_|_, 
and these have to be taken into consideration. 
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In our analysis, after including the theoretical uncertainties in the cMSSM following the 
method outlined in we set the 95% C.L. range for the branching ratio to be 



2.05 X lO"'^ < BR(Brf Xsj) < 5.05 x 10" 



(24) 



It turns out that the — )• Xs7 constraint is also extrernely sensitive to the sign of fi. For 
/i < 0, it eliminates a large part of the parameter space [43], while for fi > the constraint 



is comparatively weaker. Neither of these constraints is as strong as those arising from the 
muon (g — 2)/2; however, they are complementary to it. 

The leptonic decay Bg — )• fj,^i^i^: Within the SM, the fully leptonic decay Bg — )• is 
chirally suppressed; the SM prediction is 



BR(B, ^ 



jSM 



(3.19 ±0.35) X 10 



-9 



(25) 



The uncertainty in the BR comes principally from the decay constant fs^ = 238.8 it 9.5 MeV 



18( 1 and from the CKM element \Vts\ = 0.041 zt 0.001 [35|]. The current experimental upper 



bound by the CDF Collaboration is 



44| 



BR(B3 ^ /U+/i")cDF < 4.3 X 10"^ (95% C.L.) . 



(26) 



After including the theoretical uncertainties, we get the 95% C.L. upper limit 



BR(B, 



l-L-^H') < 4.8 X 10" 



(27) 



Inclusion of charged Higgs bosons, whose left- and right-chiral couplings depend on cot (5 
and tan/3, respectively, has a direct impact on the BR for Bs — s- which gets enhanced 

considerably above the SM prediction, and easily saturates the upper bound, especially for 
ow values of M_(_ and large tan/3. Indeed, for large tan/3, the cMSSM contribution is known 



45( 1 to scale as tan^ /3/(Ml — MH,)^. Thus, this process also constrains MFV models with a 



charged Higgs boson. It turns out that for the cMSSM, these constraints are not more severe 
than the combination of all other constraints; however, we shall demonstrate later that they 
do have an impact if the assumptions of the cMSSM are relaxed. 

In this section, we have listed the major constraints, apart from the new data on B^ — t- t+z^,-, 
on the parameter space of SUSY models, of which the cMSSM will be showcased in the following 
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section. It may be noted in passing that this hst is not fixed for ah time, as there are several other 
low-energy processes and direct search bounds which also constrain the SUSY parameter space. 
Current data on these rule out patches of the parameter space which are subsumed in the disallowed 
regions arising from the constraints which we have listed above. However, it is entirely possible 
that a future measurement - including some LHC searches - could rule out wider patches of the 
SUSY parameter space, and then the relevant processes would have to be taken into consideration. 
With this caveat, we now turn to the explicit constraints on the cMSSM parameter space, and the 
impact of the B'^ — t- t^z^t- measurement on this analysis. 



IV. CONSTRAINING THE CMSSM 



As mentioned in the Introduction, the parameter space of the cMSSM has four unknowns, viz. 
mo, ?7ii/2) ^0) tan/?, that can take real values. The sign of the fi parameter is also undetermined, 
but as indicated in the previous section, the muon (g — 2)/2 constraint disfavors < 0. We 
therefore restrict ourselves to ^ > and, hence, consider a simply connected parameter space of 
four dimensions. 

The theoretical ranges of the parameters mo, mi/2i and ^o are, in principle, completely undeter- 
mined, but the region of interest is clearly that which would lead to sparticle masses kinematically 
accessible to current accelerators such as the LHC. Keeping this in mind, we scan the ranges 

< mo < 2 TeV , < < 1 TeV , -2 TeV < ^o < 2 TeV . (28) 

The range of the remaining parameter tan /? is determined mainly by its impact on the scalar Higgs 
sector of the cMSSM, where, indeed, it arises. For very low values of tan/3 (~ 1), one tends to 
predict the lightest Higgs boson to have a small mass, which is already ruled out by the LEP 
constraints. On the other hand, if tan/3 > m^/m;,, the couplings of the charged Higgs boson to a 
tb pair begin to enter the nonperturbative regime. We have chosen, therefore, the reasonable range 

4 < tan /3 < 50 . (29) 

Using these parameters, we perform a numerical scan over the cMSSM parameter space, using 
(a) SuSpect [40] to generate the mass spectrum (this also takes care of the theoretical and direct 
search constraints), (6) SuperIso to calculate the variables listed as indirect constraints in the 
previous section, and, finally, (c) micrOMEGAs 



48( to calculate the dark-matter relic density. 



All of these are state-of-the-art software in the public domain, guaranteed to include higher-order 
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corrections as available at the moment, which have been tested in multifarious applications, as the 
literature testifies. 
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FIG. 2. The 95% C.L. constraints on the mo-mi/2 plane, for Aq — and tan/3 = 10 (left panel) and 
tan /? — 20 (right panel). Horizontal (vertical) hatching indicates regions ruled out by the measurement of 
the muon {g — 2) /2 [BR{B^ T^Vr)\ and the cross-hatched region represents their overlap. The constraint 
from B — !■ is subsumed in that from the lower bound on the Higgs mass from direct searches, and hence 
is invisible in these plots. We take /i > in this and all subsequent plots. 



In Fig. O we present the constraints on the cMSSM parameter space in the mo-mi/2 plane, 
where the tension between various measurements appears quite clearly. For this figure, we have 
set ^0 = 0, which is an assumption commonly made for simplicity. We take ;U > as required by 
the muon i^g — 2)/2 constraint. The panel on the left (right) corresponds to tan/3 = 10 (20). The 
dark areas are ruled out by the theoretical and direct search constraints explained in the previous 
section, with the approximate areas highlighted in white letteringj. Focusing on the left panel, it 
is clear that the muon (^g — 2)/2 constraint (indicated by horizontal hatching) is very stringent, 
ruling out almost all the region considered, and allowing only a small patch with low values of 
mo and mi/2- This patch, however, is disallowed by the new measurement of B'^ — t- t^Vt^ as can 
be seen from the vertical hatching. The overlap between the disallowed regions is indicated by 
cross-hatching. It is quite clear, therefore, that for tan /? = 10, there is no region in the parameter 
space shown that is at once consistent (to 95% C.L.) with both the anomalous muon magnetic 



^ Following common practice, we do not delineate separate patches in the firmly disallowed (dark) region in detail, 
as that would not be germane to the present discussion. 
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moment and the — )• r+z^,- branching ratio. We have checked that even if we take tan/3 down 
to values as small as tan/3 = 4, the two measurements, taken together with the firm constraints, 
do not allow for a simultaneously allowed parameter space. At higher value of tan /3, the situation 
is even worse. This is apparent from the right panel, where our results are plotted for tan/3 = 20. 
Here it is true that a larger region is permitted by the muon [g — 2)/2 measurement, but the region 
disallowed by the — )• r^Vr branching ratio is also much larger and covers the entire region 
allowed by (^f — 2) /2. The region disallowed by — )• t~^Vt- grows for larger value of tan /3, as may 
be guessed from Fig. and for values of tan /3 ~ 50, it would cover the whole of the parameter 
space shown in the panels of Fig. [2l Thus, for Aq = 0, one may say that these two measurements 
alone are enough to ensure that the full mSUGRA parameter space is strongly disfavored. 

It may be noted in passing that among the other constraints from the low-energy data, the 
patch disallowed by B^ — t- Xg^ is subsumed in that from the Higgs-mass bound, and, likewise, 
the patch inconsistent with Bg — )• /u^/i^ is overlaid by the dark region corresponding to the firm 
constraints. We have not, therefore, shown these disallowed regions in Fig. [2l 

The above result, disappointing as it may appear, is by no means the end of the road for the 
cMSSM, for it has been obtained only on the slice of parameter space for which Aq = 0. The 
situation changes when we permit Aq to vary. This affects the running of the charged Higgs boson 
mass, and we find that for large negative values of Aq, for a given tan/3, the Higgs-mass M_|_ is 
driven to larger values than what one would obtain by setting Aq = 0. In the context of Fig. Wib), 
this effect then pushes the model horizontally in the M_(_-tan /3 plane, eventually penetrating into 
the "allowed" region. The B^ — t- r^Vr constraint, therefore, can be quite considerably weakened 
by choosing large negative values of Aq. We do not expect such a significant change in the muon 
{g — 2)/2 constraint, but some relaxation is not unreasonable to expect when Aq is varied over a 
wide range. Large negative values of also tend to increase the mass of the lightest Higgs boson 
, thereby relaxing somewhat the LEP bounds arising from Higgs boson mass considerations 49l |. 



Accordingly, we repeat our analysis of the constraints on the raQ-raii2 plane, keeping floating 
between —2 TeV and +2 TeV. In our analysis, a point in the mo-?Tii/2 plane, for a given tan/3, is 
taken to be allowed at 95% C.L. by a given constraint if we can find any value of Aq, lying in the 
range —2 TeV < ^ 2 TeV, for which the given constraint is satisfied. Our results are exhibited in 
Fig. [21 which follows the notations and conventions of Fig. [2] closely. In addition to the constraints 
shown therein, the one from B^ — )■ X^'^ now makes its appearance as a small region hatched with 
slanting lines, indicating that this constraint is now stronger than the weakened Higgs-mass bound 
from collider machines. However, it is not strong enough to rule out any portion which is not 
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already disallowed by the other constraints. 
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FIG. 3. The same as in Fig. [21 except that the trilinear coupling Aq is kept floating over the range 
(—2, +2) TeV, and regions are considered disallowed only if they remain disallowed for all values of Aq 
in the given range. This weakens the constraints enough for a small allowed region to appear in the left 
panel (tan/3 = 10). The yellow/light gray region is simultaneously consistent with all constraints at 95% 
C.L. 



Even a cursory examination of the dark and hatched regions in Fig. [3] will indicate that, while 
the qualitative features of the regions disfavored by the muon (g — 2)/2 and the — t- t'^Vj- 
measurements stay the same, somewhat larger areas in the plane are "allowed" by each constraint 
individually. This, by itself, is not surprising, but it has the exciting consequence that now the 
left panel (tan /3 = 10) exhibits a small patch, roughly triangular in shape, which satisfies each of 
the constraints individually for some value of and moreover, there is a subset of this region 
where all the constraints are satisfied simultaneously for the same value of ^o- This subregion, 
which represents the actual parameter space consistent with all the measurements individually at 
95% C.L., is denoted by yellow/light gray shading. It is on this "allowed" region that we focus our 
interest in the subsequent discussion. 

If we glance at the right panel of Fig. [3l where tan/3 = 20, we see that there is no allowed 
region at all, the disallowed regions showing substantial overlap and covering the whole of the 
plot area. Once again, we surmise that high values of tan/3 are disfavored, whatever value of Aq is 
chosen, and that the allowed region in the cMSSM parameter space must lie in the neighborhood of 
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tan/3 = 10. We have already mentioned that consistency with the — )• t'^i't constraint requires 
large negative values of the Aq parameter, which can drive Af+ to higher values even for the low 
tan/3 (w 10). In order to see this, we plot, in the left panel of Fig. HI the same constraints in 
the plane of Aq and "1-1/21 keeping mo fixed at the value tjiq = 150 GeV, for tan/3 = 10. This 
particular value of mo has been chosen since in the left panel of Fig. [3l it lies more-or-less near 
the center of the allowed triangle (yellow/light gray shading) and is roughly the value for which 
the maximum range of mi/2 appears to be allowed. The dimensions of this allowed triangle also 
encourage us, in Fig. [H to "zoom in" on the range mx/2 = 300-500 GeV, outside which we get 
disallowed regions. However, ^0 is varied between —1.5 TeV and +1.5 TeV, to adequately cover 
the whole range allowed by the firm constraints, as is apparent from the left panel of Fig. [H 
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FIG. 4. Further constraints on the cMSSM parameter space. The left panel shows constraints on the mi/2- 
Aq plane, for mg = 150 GeV and tan/? = 10. The right panel shows, similarly, constraints on the tan/3-Ao 
plane for too = 150 GeV and mi/2 — 400 GeV. Notations and conventions are the same as in Figs. [2] and [H 
The dotted (red) line represents the dark-mattcr-compatible region, and the black dot superposed on it is a 
benchmark point chosen for LHC studies. 



It is immediately apparent from the left panel of Fig. U] that about half of the region with 
positive values of Aq is ruled out by the firm constraints, and the remaining half by the — )• t'^i't 
measurement. The latter has a severe impact on the Aq < region as well, essentially forcing us 
to consider large negative values of Aq for small values of mi/2- Including the muon (g — 2)/2 
constraint, which disfavors large values of mi/2) then clinches the issue, permitting only another 
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small wedge-shaped (yellow/light gray) region allowed by all the constraints. The maximum range 
of mi/2 permitted by all the constraints is around 325-425 GeV, which matches tolerably well with 
the vertical limits of the allowed triangle in Fig. [31 as should be the case. The value Aq, on the 
other hand, is quite strictly restricted to the approximate range —625 GeV to —1.4 TeV. 

Of course, the above results are only for a fixed value tan (3 = 10. Though we have already seen 
that jumping to a much larger value tan /3 = 20 does not lead to any allowed region, it is interesting 
to 'zoom in' to the tan ^-Aq plane and see the impact of all these constraints there. This is shown 
in the right panel of Fig. U where we set mg = 150 GeV as before, and mi/2 = 400 GeV, which 
lies close to its value on the left panel for which the allowed range of Aq is maximum. Once again, 
the combined constraints predict a small allowed region, with a maximum range of tan/3 lying 
roughly between 8 and 12. If we now refer to Fig. ^b), this means that the charged Higgs boson 
is predicted to have a mass Af+ > 600 GeV. 

Combining all these results, therefore, we obtain a roughly polyhedral allowed region in the 
four-dimensional parameter space, which is enclosed in a rather small hypercube with sides ap- 
proximately at 

100 GeV < mo < 225 GeV , 
375 GeV<mi/2< 425 GeV, 
-1.4 TeV < Ao < -625 GeV , 

8<tan/3<12. (30) 

The volume of the actual allowed region is considerably smaller than that of the hypercube, given 
that the two-dimensional projections shown in the previous two figures are roughly triangular in 
shape. Compared to the large regions considered in traditional work on the cMSSM, this constitutes 
a rather specific region of parameter space, and encourages us to make specific predictions based on 
this model. One can easily argue that the qualitative features of the mass spectrum and couplings 
will not undergo dramatic changes from one end to the other of so small a box as this one, unless, 
indeed, it encloses some point(s) of instability. This is unlikely, for none of the many studies of the 
cMSSM parameter space have ever shown such a possibility. 

The very first prediction one would naturally demand from a specific point or region in the 
cMSSM parameter space is whether this can adequately explain the dark-matter content of the 
Universe relic density of LSP's. The CMBR data indicate a relic density Qh^ = 0.1123 ±0.007 
at 95% C.L. [13]. In general, SUSY models with a low- lying mass spectrum, such as the one 
in question, tend to predict too large a density of LSP's unless these are coannihilated by some 
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reaction with a substantial crosssection. This leads to a restriction on the cMSSM parameter 
space, which, given the accuracy of the CMBR data, confines us, more or less, to a line passing 
through the four-dimensional parameter space. The dark-matter requirement is known to favor 
large negative values of 13 ; it is rather gratifying to see that this line passes right through 
the allowed region in the parameter space discovered in this work - which seems to indicate that 
a SUSY explanation of dark-matter may indeed be the correct one. The line consistent with the 
dark-matter requiremenljll is shown by red dots on both panels in Fig. U] and may be seen to pass 
clearly through the allowed region, favoring a narrow range of Aq around —1.25 TeV and tan/3 in 
the range 9 — 11. This region is rather close to the forbidden stau LSP region, as is apparent in 
both panels of Fig. [H In the allowed region, the lighter stau ri is rendered very light due to the 
presence of a large negative nirAr in the off-diagonal terms in the stau mixing matrix; however it is 
marginally heavier than the neutralino x? LSP. This permits stau coannihilation with neutralinos, 
and reduces the relic density so that it is within the observed range. 

We note that there does not seem to be any a priori reason for the region allowed by the 
low-energy constraints to match with the dark matter-compatible region, since the low-energy con- 
straints come from processes quite different from those that control the relic density. Nevertheless, 
the fact that the two regions do show some overlap encourages us to argue that we are now con- 
verging on the correct region in the hitherto-unknown parameter space. We may, therefore, make 
bold as to venture some predictions regarding the collider signals for this range of parameter space, 
especially in the context of the LHC. 

In order to make a clear prediction about the LHC signals, we choose the following benchmark 
point in the cMSSM parameter space: 

mo = 150 GeV , mi/2 = 400 GeV , = -1250 GeV , tan/3 = 10 , ^ > . (31) 

Not only does this lie inside the hypercube marked out in Eq. ()30p . but it lies well within the 
allowed region, and right on the line corresponding to the dark-matter requirement. In Fig. [H this 
benchmark point is indicated by a small black circle in both panels [on top of the dotted (red) line 
labeled "dark matter"]. The mass spectrum and signals expected for this "golden point" will be 
typical of the entire allowed region, which is, after all, rather small. At this benchmark point, we 
get the central values of the observables to be BR{B Xsj) = 2.64 x IQ-"^, R!^f = 0.993, and 
= 13.0 X 10^^'^. Clearly, all of these are consistent with the measurements to within 2a, though 
R^J" only barely survives the 2a bound. 

^ We do not go so far as to call it a constraint, though this is not unheard of in the literature. 
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Let us first discuss the cMSSM mass spectrum expected with this benchmark point. We cal- 
culatqj the mass spectrum and the branching ratios using the code SUSY-HIT 5l|] and taking 
rrit = 173.1 GeV. The Hghtest Higgs boson is predicted to have a mass around 119 GeV, which 
is consistent with current bounds, but lies precisely in the range where its detection is most prob- 
lematic because of large QCD backgrounds at the LHC. In fact, a light Higgs boson of this mass 
range must be detected through the rare decay — )• 77, which is unlikely in the 7 TeV run, and 
will require the accumulation of a fair amount of statistics even in the 14 TeV run. The heavy 
Higgs bosons, including the i/"*", will lie in the range 835-840 GeV, which is again kinematically 
inaccessible in the 7 TeV run, but may be detectable at 14 TeV. We have already shown that for 
tan /3 = 10, as taken for this benchmark point, this high value of is what allows us to evade 
the —7- T^Ur constraint. Turning now to sparticles, the LSP Xi will have a mass of 164 GeV, 
with the next-to-LSP (NLSP) being, as expected, the ti with a mass of 171 GeV. As explained 
above, the closeness in these masses permits the coannihilation of stau, so that the relic density 
is controlled. In this scenario, this stau and the lightest neutralino are the only sparticles with 
masses below that of the top quark, all other particles being heavier. The nearly degenerate lighter 
chargino xt second neutralino X2 li^ at 315 GeV, while the other sleptons and the sneutrinos 
have different masses in the 200-320 GeV range. The gluino mass, however, is as high as 934 GeV 
and the squark masses mostly populate the range 800-900 GeV, except for the 61, with mass around 
719 GeV, and a light stop ti which lies as low as 393 GeV. 

An immediate consequence of these large squark and gluino masses is that the sparticle pro- 
duction cross section at the LHC will be on the low side: at 7 TeV it will be around 0.4 pb at the 
leading order (LO), while at 14 TeV, it will have the much healthier value of 5.2 pb at LO. About 
60% of these crosssections come from squark pair production, of which roughly half is due to tit^ 
production alone. The ti will decay to a top quark and a neutralino with a BR ~ 2/3, and hence, 
a possible signal would be a top-enriched final state with large missing transverse energy (MET). 
However, the enormous background to this process must be taken into consideration when study- 
ing this signal. The other traditional signals for SUSY - cascade decays of the gluino or squarks to 
charginos and heavy neutralinos, ending up in multileptons, jets, and MET - in this case provide 
r-rich final states because of the low-lying ri. However, r's coming from the decay ti — t- r + 
will generally be too soft for detection, because of the small splitting M(ti) — M{xi) — 7 GeV. 
Final states involving other charged leptons will be suppressed. This indicates that the best option 



The masses obtained using different RG evolution algorithms differ by a few GeV, and the errors from the calcu- 
lation are difficult to quantify. Here we give the exact values obtained by SUSY-HIT. 
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to seek SUSY with this benchmark point is the final state with four or more jets and substantial 
MET, which can arise from cascade decays involving only strongly interacting sparticles and the 
invisible LSP. 

1000 

tan P = 10 

800 



600 
(3 



E 400 




|x >0 



ATLAS : 5a,7TeV, 1 fb 



Higgs mass 



500 



1000 1500 

mo[GeV] 



FIG. 5. The mo-mi/2 plane in the cMSSM, showing the allowed region (yellow/light gray) for floating Aq, 
as well as the 5cr reach of the ATLAS detector at the 7 TeV run, using the jets + MET signal. The black 
dot inside the allowed region is the golden point discussed in the text. The entire nonshaded region inside 
this plot is accessible to the LHC with ^/s = 14 TeV with a luminosity of 10 fb~^. 



In Fig. Owe show the allowed parameter space (yellow/light gray) in the mn~»T,i / 9 p lane for 
tan/3 = 10 and floating Aq (as in Fig. [3]), and also the ATLAS 5a discovery limit 64] at the 
7 TeV run with an integrated luminosity of 1 fb~^, using the four-jets + MET channel. It may 
be seen that the entire parameter space allowed by low-energy constraints at 95% C.L., including 
our golden point, lies just outside the 5a discovery limit of ATLAS. The ATLAS study, in fact, 
has shown that at neighboring points, an overall crosssection of about 1 pb is required for a 5a 
discovery using the four-jets -|- MET channel. Making a simpleminded scaling with the predicted 
crosssection of 0.4 pb at our benchmark point, one may expect a signal in this channel at the level 
of about 2a. The same conclusions can be reached using the CMS 95% exclusion plot for the 
7 TeV run [l^. Thus, by the end of 2011, as per the LHC running schedule, we may begin to see 
tantalizing hints of SUSY. If this should occur, then, in the 14 TeV run, it will be easy to see a 5a 
signal with even 1 fb^^ of data - which should collect within the first few months. Interestingly, 
some of the direct production modes for charginos, which are electroweak in nature, lie at the level 
of 5%-10% of the total crosssection. These may be difficult to detect in the 7 TeV run, but in the 
14 TeV run, they are sure to provide additional signals for SUSY. With such copious production 
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of sparticles, the LHC could indeed act as a SUSY factory, as mentioned in the Introduction. 

V. NUHM : EXPLAINING B+ ^ r+Vr 

In the previous analysis, we have seen that the combination of constraints on the cMSSM 
parameter space leads to the prediction of a small value of tan /3 and hence, according to Fig. 
the charged Higgs boson is necessarily heavy. Comparison with Fig. [I^a) readily shows that in this 
limit, the model is only just consistent with the B'^ — t- t~^Vt- constraint at 95% C.L. However, if 
we take the position that the 2a discrepancy between the SM and the experimental result should 
be explained by a positive NP contribution, then the cMSSM fails the test, for it actually tends 
to diminish the SM prediction, and barely survives exclusion in a decoupling limit. This bare 
survival, by the skin of its teeth, as it were, is the proximate cause of the stringent constraints on 
the cMSSM parameter space discussed in the previous section. 

As the cMSSM is the SUSY model with the maximum number of simplifying assumptions (and 
hence the minimum number of free parameters), it is interesting to ask if the relaxation of one 
or more of these assumptions could lead to a SUSY model which actually explains, rather than 
merely remains consistent with, the — t- t'^Vt discrepancy. Since the NP effect in i?+ — t- r^Ur 
involves the scalar sector of the cMSSM, an obvious option would be to consider a model where 
the parameters of the Higgs sector are given a greater degree of flexibility than in the highly 
constrained cMSSM. In this context, an obvious choice of model is the so-called nonuniversal Higgs- 
mass (NUHM) model, which is an extension of the cMSSM where the Higgs-mass parameters rriH^ 
and are delinked from the universal scalar mass parameter mo at the GUT scale and are 
allowed to vary freely ^]. At the electroweak scale, these two extra parameters mn^ and m/^j are 
usually traded for the Higgsino mixing parameter and the pseudoscalar Higgs boson mass m^. 
This model, therefore, has six parameters, viz. mQ,mi/2, fJ', Ma, Aq and tan/3. 

NUHM models have been studied rather extensively, and various constraints on the six- 
dimensional parameter space have been found and exhibited in the literature 
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5j]. What 



interests us here is the fact that is a free parameter in the model, and it can be easily exchanged 
for M+, to which it is related by the well-known SUSY relation 



Mi = Mi + M^, (32) 

at tree level. We can accordingly fix mo, ^711/2; ^O; etc. at whatever value is required to satisfy the 
other constraints in the cMSSM, and then claim an explanation for the B^ — )■ t'^i'j- discrepancy by 
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choosing a low M4. and a high tan /3 - this freedom being allowed by the bigger parameter space in 
the theory. However, large values of tan j3 and small values of Af+ lead to large charged Higgs boson- 
mediated contributions to the FCNC process Bg — )• thus restricting the freedom in choosing 
parameter values. Here, as explained earlier, SUSY cancellations between the charged Higgs boson- 
mediated and the gaugino-mediated contributions come to the rescue: stringent bounds can be 
evaded if the gaugino masses are somewhat low, comparable to that of the light charged Higgs 
boson H'^ . This, in turn, demands that the universal gaugino mass be somewhat small, 
compared with the other parameters, which are not so restricted. 
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FIG. 6. Constraints on the Ma-h parameter space in the NUHM model. Notations and conventions are 
the same as in the previous plots, except that now there is a significant constraint from Bg — > rather 
than Bd — > ^^7. The yellow/ light gray region is allowed by all the constraints, and the black dot inside it 
is a benchmark point chosen for LHC studies. 



While an exhaustive study of the NUHM parameter space vis-a-vis the present set of constraints 
would require a separate work in itself, it is interesting to see if the NUHM model can at all provide 
regions in parameter space which are consistent with all the constraints, and can simultaneously 
provide a NP explanation of the — t- r^Vr discrepancy. To illustrate that this is, in fact, possible, 
we show in Fig. [6] the regions allowed by the different constraints in the plane, keeping all 

the other parameters fixed at 



mo = 1.2 TeV , 



"1-1/2 



180 GeV 



A. 



1.2 TeV 



tan /3 = 50 . 



(33) 



A glance at the figure will reveal that here, as in the cMSSM, there is a complementarity between 
the B^ —7- T^v-r constraint and the [g — 2)/2 constraint, the former tending to rule out larger 
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values of Ma and the latter tending to rule out larger values of /i, as a result of which only a 
small rectangular patch in the rriA-n plane is allowed by both constraints taken together. A large 
portion of this remaining patch is again disallowed by the Bg — )■ constraint, leaving a roughly 

sickle-shaped yellow/light gray region. In this region, Ma remains in the approximate range 100- 
150 GeV; i.e. M_|_ lies roughly in the range 125-170 GeV, according to Eq. ()32|) . Figure [1] then 
tells us that this is not only consistent with the experimental data, but it is precisely the range 
for which the NP explanation saturates the gap between the SM prediction and the experimental 
central value. 

As before, to be precise about the LHC signals, we choose a benchmark point, which has the 
fixed parameter choices of Eq. ([331) as well as 

Ma = 145 GeV , // = 400 GeV , (34) 

which is indicated in Fig. [6] by a small black dot in the middle of the allowed (yellow/light gray) 
patch. The central values of the observables at this point are BR(B — t- Xg^) = 3.50 x 10~^, R^J' = 
1.24,BR(Ss /u+/i") = 3.22 x 10"S,a^ = 12.8 x 10"^°, ah of which are weh within the 2a range 
of the respective measurements. We note that the relic density of LSP's at this point is not enough 
to saturate the CMBR requirements, which means that this model is not ruled out by the latter, 
but is not a solution to that problem either. 

The major features of the mass spectrum at this benchmark point are as follows: the lightest 
Higgs boson lies just beyond the LEP disallowed region, at 112 GeV, and, as in the cMSSM, this 
is a difficult mass range to search for the lightest Higgs boson. We would have to wait for enough 
statistics to accumulate at the 14 TeV run to see this Higgs boson in the 77 channel. The 
and lie at 145 GeV and may just be detectable through their decays to WW* modes, while 
the charged Higgs boson lies at 170 GeV, where it will decay to t^Ut- These may also be 
detectable fairly early in the 14 TeV run. The LSP, as before, is the lightest neutralino xi with a 
mass of 71 GeV, which is permitted by the LEP direct search bound as applied to the NUHM 351 ]. 



The Xi s-iicl X2 lie at around 130 GeV, while the other gauginos are heavier than 400 GeV. The 
sleptons and squarks in this model are very heavy, lying in the range 700 GeV to 1.2 TeV, but the 
gluino 'g is comparatively light, having a mass of 511 GeV. 

As a consequence of the low-lying gaugino states contrasted with heavy sfermions, the dominant 
sparticle production channels in this model turn out to be to chargino pairs xtxi (~ 50%) and 
chargino- neutralino pairs xf X2 (~ 25%), with gluino pairs Ijg bringing up the rear (~ 20%). The 
total crosssection in this model would be 4.4 pb at 7 TeV and 27.3 pb at 14 TeV, i.e. much larger 
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than the earher case of the cMSSM. The gluino production channel can give rise to the same jets 
+ MET signal as before, as each gluino will undergo three-body decays through virtual squarks. 
The production crosssection for 'gg pairs at 7 TeV is around 0.8 pb, which indicates that the jets + 
MET signal may actually be observable in the 7 TeV run at the 3a-Aa level when 1 fb"-*^ of data 
have been collected. At 14 TeV, of course, a few hundred pb^^ of data would be enough to obtain 
a 5(7 signal in this channel. Turning now to the chargino production modes, the rate of production 
of XiXi indicates a crosssection for a dilepton + MET signal around 90 fb, which may not be 
discernible above the background, especially as the mass splitting M{xi) — -^(Xi) is rather small. 
However, the X1X2 channel could lead to hadronically quiet trilepton + MET signals at the level 
of 74 fb, which have smaller SM backgrounds and hence could probably be seen as more data are 
collected in the 7 TeV run, and would be a sure-shot option at the 14 TeV run. 

Before concluding this section, we reiterate that NUHM models which explain the — t- r^i^r 
discrepancy and, at the same time, remain consistent with the data on Bs — )■ fi^fi^, will generically 
come with light gauginos, and lead to collider signals somewhat similar to those discussed above. 
However, what we have studied is just one portion of the NUHM parameter space, inasmuch as 
we fixed mo and Aq to very large values. A more comprehensive scan over the NUHM parameter 
space might reveal more patches consistent with all the constraints, and some of these may lead 
to collider signals which are different from those discussed in the context of our benchmark point. 
The detailed exploration of the NUHM parameter space in this context calls for a separate study. 

VI. CONCLUDING REMARKS 

With the commissioning of the LHC, the search for new physics beyond the Standard Model 
has assumed paramount importance in particle physics at the high scale. However, low-energy 
observables from flavor physics, like those from the decays of K, D, or B mesons, can offer indirect 
constraints on high scale physics. Indeed, with the high statistics available at the B factories 
BaBar and Belle, the freedom available for new physics has been substantially constrained. Most 
of the low- energy measurements have been consistent with the SM, and hence allow only a little 
leeway for NP. On the other hand, it is seen that the handful of measurements that indicate & ^ 2a 
deviation from the SM also restrain the NP parameters from taking arbitrary values. 

In this paper, we have shown that the combined effect of both these kinds of low-energy measure- 
ments - those consistent with the SM (e.g. the branching ratios of B^ Xg^ and Bg IJ-'^IJ-~) as 
well as those showing deviations from the SM (e.g. the anomalous magnetic moment of the muon. 
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and the branching ratio of the — )• T'^i'r) ~ results not only in indicating new physics, but also in 
pinpointing the relevant new physics parameters. In particular, we have pointed out that the latest 
measurement of — )• t^Ut branching ratio has a large impact on a large class of NP models, 
especially those which include a charged Higgs boson H^. In fact, the decay — )• t^Ut, by itself, 
can constrain most of the models with minimal flavor violation that involve an . This is because 
the latest measurement gives a branching ratio ~ 2a more than the SM prediction. If this discrep- 
ancy is to be explained by a MFV model, one needs very light charged Higgs bosons < 200 
GeV) and large tan/3 (> 20). On the other hand, a heavy charged Higgs boson (M-(_ > 300 GeV) 
and a small tan /3 can be barely consistent with the data to within 2a, but cannot be considered 
an explanation for the gap between theory and experiment. This is a general result that can be 
applied to any member of the MFV models, and we choose to apply it to the constrained MSSM, 
which is motivated by mSUGRA and is one of the most predictive SUSY models. 

In cMSSM models, the charged Higgs boson is typically heavy, so that only the low tan /3 region 
survives the B^ — t- r"*"!/,- measurement. When combined with the anomalous magnetic moment 
of the muon, the fate of even this region is in jeopardy: indeed, for a vanishing universal trilinear 
coupling Aq, there is no region in the cMSSM parameter space that is consistent with both these 
measurements to 95% C.L.. The situation can only be salvaged with a large and negative and 
that too for an extremely small region in the mQ~mi/2 plane. The combined low-energy data thus 
pinpoint us to a very specific location (the golden point) in the five-dimensional parameter space 
of cMSSM: /i > 0, Aq ^ -1.25 TeV, tan/3 10 ,mo ~ 150 GeV,mi/2 ^ 400 GeV, with a spread 
not more than that given in Eq. (|30p . It is remarkable that for part of this specific region, including 
the golden point, the mass and coupling of the LSP are exactly such that it can account for all the 
dark-matter in the Universe. This may either be a coincidence, or an indication that we are on the 
right track in our quest. 

If we indeed are on the right track, and the golden point of the cMSSM is actually the NP that 
we have all been looking for, then we may not have to wait too long for its discovery. Since the 
values of rn-o and ^f this point are rather small, at the LHC, one expects a weak 2a signal 

in the jets -|- MET channel even in the 7 TeV run with 1 fb^^ of integrated luminosity, and a 5a 
discovery early in the 14 TeV run with just 1 fb~^ of data. 

While the above suggestive coincidence is quite appealing, and the prospects of the detection 
of SUSY during the early parts of the 14 TeV run quite enticing, even at the golden point the 
model barely survives the 95% limit bounds and does not offer any help at all in explaining the 
B^ —7- T^i'r data - the absence of a light makes it impossible for the cMSSM to do so. We 
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therefore explore a related but less constrained model, the NUHM, where the charged Higgs boson 
mass can be considered to be a free parameter. Here the presence of an extra parameter works 
wonders for explaining the low-energy data, covering the entire experimentally allowed region, 
including the central value. This model can also lead to rather spectacular trilepton + MET 

signals at the LHC, which may become detectable even towards the end of the 7 TeV run. 



While our work tends to indicate a rather specific region of parameter space and specific sig- 
nals at the LHC, especially for the cMSSM, there are some caveats which need to be taken into 
consideration, even apart from theoretical issues in the construction of the cMSSM. The first is 
the issue of experimental errors on the low-energy measurements, which we have taken at the 2a 
level. If these are given more latitude (e.g. taken at the 3a level) the constraints from low-energy 
processes would be considerably relaxed. In particular, the S+ t'^Vj- measurement would still 
allow wide regions in the cMSSM parameter space. However, it would still disfavor very large 
values of tan /5 ^ 50. A more serious point is the asymptotic behavior B^^^ — t- 1 in the large M_|_ 
limit, as compared to the 2a bound B^^^ > 0.99. The strong constraint on NP comes because 
one must squeeze the contribution of the charged Higgs bosons into the narrow region 0.99 — 1.00. 
A small downward revision in the lower bound on -R^^ could allow large tan/3 values even for 
large M+. Such deviations can come from a variety of sources, such as higher-order corrections, 
a slightly-changed value of /b or \Vujj\, or a revised experimental result. On the other hand, a 
small upward revision of the allowed R^^_^ band could rule out the entire gamut of MFV models 
with M_|_ > 200 GeV. In particular, even the small leeway allowed for the cMSSM would then be 
closed. Wc note, therefore, that the bounds and predictions presented in this paper are specific to 
the experimental limits as they stand at present. 



We have not made a very detailed study of the LHC signals, confining ourselves to generalities, 
because it is somewhat premature, at this stage, to make very definite predictions in this regard. 
Nevertheless, our work has highlighted the fact that if indeed we are to accept the cMSSM at face 
value, as most LHC studies do, then we should take the cMSSM in its entirety, i.e. all constraints 
from all sectors, including the low-energy sector. The next year and the years after it will be the 
most crucial in determining if our analysis, in fact, is on the right track. 
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